Photogalvanic current in a double quantum well 
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We study the in-plane stationary current caused by phototransitions between the states of a double 
quantum well. The electric polarization of light has both vertical and in-plane components. The 
stationary current originates from the periodic vibration of electrons between two non-equivalent 
quantum wells caused by the normal component of the alternating electric field with simulteneous in- 
plane acceleration/deceleration by the in-plane component of electric field. The quantum mechanism 
of the stationary current is conditioned by in-plane transition asymmetry which appears due to the 
indirect phototransitions with the participation of impurity scattering. The photocurrent has a 
resonant character corresponding to the equality of the photon energy to the distance between 
subbands. It is found that the current appears as a response to the linear-polarized light. 

PACS numbers: 73.50.Pz,78.67.-n,72.40.-|-w, 73.21. Fg 



I. INTRODUCTION 

Since the first studies on the photogalvanic effect 
(PGE) at the end of the 70th, a wide literature devoted 
to this subject has appeared [J-Q, see also reviews [^-Q- 
The activity in this field continues up to now (see, for 
example, |l(]| - [l8| ). There are different variants of PGE 
in confined systems: the stationary in-plane photocur- 
rent in classical [l^ and quantum feo, 21] films, and the 
current along solid-state surface |22l424{ . This photocur- 
rent exists even if crystal asymmetry is negligible, but 
the quantum well is oriented (directions across the well 
are not equivalent). The current along the surface occurs 
if the electric field of the light has both in- and out-plane 
components. 

The phenomenology of PGE is determined by the re- 
lation for the current density 

j = a,((E - n(nE))(nE* + c.c) + iaa[n[EE*]], (1) 

where n is the normal to the quantum well. Real con- 
stants Us and Ua describe linear and circular photogal- 
vanic effects, correspondingly. The origin of this current 
can be understood if to consider the out-of-plane com- 
ponent of electric field as modulating the quantum well 
conductivity with a simultaneous driving of electrons by 
the field in-plane . 

In a quantum well the vertical component of the elec- 
tric field of light can cause the transitions between dif- 
ferent quantum subbands. In the presence of scattering 
this gives birth to the effective pumping of the in-plane 
momentum to the electronic system. The light plays the 
role of the energy and non-equilibrium source, while the 
scatterers produce electrons in-plane acceleration. The 
situation is, in a certain sense, similar to the motion of a 
car where the friction forces the car to move. 

The purpose of the present article is to study the mech- 
anism of PGE in a double quantum well. This system 
looks perspective because the structure of the levels of a 
double quantum well permits easy tuning of the distance 
between subbands to the frequency of the external field. 



The effect under consideration is illustrated in Fig.l. 
We consider intersubband transitions of electrons in a 
system with the quadratic energy spectrum. An electron 
goes between two states en(p) and e„'(p') due to the si- 
multaneous action of electric field and scattering. These 
states originate from mixing the states of different indi- 
vidual quantum wells. The in-plane current appears due 
to the change of electron in-plane momentum. To " mem- 
orize" electric field in- and out-plane components, the 
transition probability should contain their product. For 
non-conservation of the electron momentum the scatter- 
ing should be taken into account. This transition proba- 
bility arises in the second order of the perturbation the- 
ory. The amplitude of transitions has a resonance on 
an intermediate state. The subbands of the quantum 
well are equidistant, that gives rise to the absence of the 
resonance smearing due to the difference in electron mo- 
menta. The result of excitation is the pumping of the 
momentum to the electron subsystem and the in-plane 
current. The paper is organized as follows. First, we will 
discuss a simple classical model of the effect based on a 
parabolic well. Then, we will find the transition proba- 
bility in a classical electric field. After that the current 
will be found using many-band kinetic equation. 



II. SIMPLE CLASSICAL MODEL 

To explain the physical origin of the effect we con- 
sider a simple classical model instead of a 2D system: an 
electron in an oscillatory well in z-direction with confin- 
ing potential mwQZ^/2 affected by the alternating electric 
field with x and z components E(i) — Re(Ee~''^*). The 
classical Newton equation for an electron reads 



r + 7r = eE/m, 



(2) 



where we introduced the liquid friction coefficient 7 = 
7o + 71-2^- The dependence of the friction on z takes into 
account the assumed weak asymmetry (70 3> 7iz) of the 
well in z-direction. 
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FIG. 1; (Color online) a) The sketch of the proposal experi- 
mental setup. The electric field of light E(f) is tilted in (x,z) 
plane. The stationary current is directed along the x-axis. b) 
The sketch of the band structure. Quantum wells are cen- 
tered in planes z — ±d/2. The carriers are provided by the 
5-layer of donors in plane zq. c) The transition amplitude in- 
cludes vertical transition caused by light between ± subbands 
and impurity scattering which does not conserve the in-plane 
momentum. 



circular polarization. The origin of this behavior is ex- 
plained by the character of the electron motion in the 
zero approximation. Indeed, if 71 — 0, for linear polar- 
ization, the electron rotates in the exact resonance and 
vibrates along a straight line out of resonance. For cir- 
cular polarization the behavior is opposite. 

Liquid friction force — 71- does not affect the direction 
of vibrating motion; therefore it does not produce a drift. 
At the same time, due to 71, a rotating particle differently 
brakes at the opposite (upper and lower) sides of the 
circle that produces a translational displacement, and as 
a result, the mean drift. In the case of circular-polarized 
light, the direction of the motion depends on the sign 
of polarization and the sign of resonance detuning. The 
value of the drift velocity near resonance does not depend 
on the friction strength, but it depends on ratio 71/7. 

Fig. 2 illustrates this reasoning by the exact solution of 
the Newton equation. The photogalvanic effect in this 




FIG. 2: (Color online) Solution of the Newton equation 
for linear (c,d) and circular (a,b) polarized electric field for 
initial conditions r(0) = 0, r(0) = and parameters 70 = 
0.2, 71 = 0.02, eE^/m = 0.3, eE^/m = 0.1, ujq = 1; for plot 
a.) uj = 0.9, for plot b) a; = 1.1, for plots a) and h) uj = 1. 
For circular polarization the sign of detuning determines the 
direction of the steady-state drift. 



The forced solution of the Newton equation is found 
by expanding in powers of 71: 



zq — Re 



r = ro + ri 



m(— (jj'^ — i^QOj) 



model has a purely classical nature. In particular, the 
circular PGE does not need the spin pumping as in spin- 
related circular PGE. At the same time, the classical and 
quantum photogalvanic effects have different properties. 
The photogalvanic effect on intersubband transitions of 
double quantum well Considered below has a resonant 
character like the classical PGE discussed here. The dif- 
ference is the absence of circular photogalvanic effect for 
transitions in a double quantum well. 



xq = Re 



m{—uj'^ — ijo'^) 



III. TRANSITIONS BETWEEN SUBBANDS OF 
DOUBLE WELL 



(ii) 



Im 



e:e.. 



27om2 (aj2 — i7oa;)(a;2 — Wg -I- «7ow) ' 



(ii) = 0. 

_ (3) 
Here (...) denotes the time averaging. 

Let damping 70 be also small. Then the mean velocity 
has a resonance at uj — ujo- The frequency behavior near 
this point depends on the kind of electromagnetic field 
polarization: delta-like peak for the linear polarization 
and antisymmetric Fano-like resonance cx 1 / {uj — ujq) for 



We study electrons with a parabolic isotropic energy 
spectrum in a double quantum well (see Fig.l). The am- 
plitude of transition between wells is weak, but compara- 
ble to the separation of energies of individual wells. The 
states with the in-plane electron momentum p and the 
subband number n = ± |n,p >= Xn{z)exp{ipp)/^/S, 
{S is the system area, we set ?i = 1 throughout this sec- 
tion besides the final expression) have energies e„_p = 

/2m -t- e„. In this case, the subbands are parallel. 



— This circumstance plays an im- 



portant role in the further consideration, providing the 
resonance of optical frequency with a distance between 
subbands for electrons with arbitrary momenta. The 
overlapping of wave functions Xn{z) is supposed to be 
weak and intersubband distance £+ — £_= A (A > 0) is 
small as compared to the Fermi energy. The scatterers 
(donors) are distributed in a delta-layer at z = zq > 0. 
The well widths and the distance d between them are 
assumed to be small as compared to Zq. 

Assuming that the mean free time is large as compared 
to the distance between the levels of quantum wells (and 
also the Fermi energy) one can treat n and p as good 
quantum numbers and describe the problem within the 
kinetic equation for distribution functions fn,p- In such 
an equation, external classical alternating electric field 
E(t) = Re(Eoe^*"*) causes the transition between un- 
pertubed states and determines the generation term in 
the kinetic equation. The interaction with charged im- 
purities provides the mechanism of electron scattering. 
The kinetic equation reads 



E 

n',p' 



n,p;n 



,p v-' n ,p 



_ f(i)- 



Gn.p — 0, 



(4) 



where the generation Gn^p is given by 

G„.p=E<,P;«%P'(/<P'-/^°^)- (5) 

n' ,p' 

Here W™^.^, p, is the impurity transition probability, 
^np-n' p' transition probability due to the com- 

bined action of electromagnetic field and impurities, f^p 

is the equilibrium distribution function and fn,p is the 
first correction to the distribution function in the ex- 
ternal electromagnetic field. Using the classical kinetic 
equation means neglecting the off-diagonal elements of 
the density matrix that is valid if the collision broaden- 
ing of subbands is less than the distance between them. 
The perturbation includes the Hamiltonian of the inter- 
action with electromagnetic field H^'^ and the potential 
energy of the electron interacting with impurities V. The 
first is 



HP'' ^ -Re(Ae" 
c ^ 



+ h.c.), (6) 



where Re(Ae~'"*) is the vector potential of electromag- 
netic field with frequency uj, v = (vll,w^) is the veloc- 
ity operator. The complex amplitude of electric field 
is E = iujA/c. Thus, the operator U = e(Ev)/iuj. 
Note that we suppose the electric field to be homoge- 
neous. The diagonal elements of in-plane components 
of the the velocity operator v|[ p.^^, p, — Vp(5nn'(5p_p', 
p/m. The normal component has matrix 



elements 
reads 



n,p;n ,p 



n,n''^P,P'- 



The impurity potential 



y(r) = ^^^(r-r,), 



(7) 



where the sum runs over all the impurities situated in 
points Ti) with individual potentials u(r — r^). 

The appearance of the photogalvanic current requires 
non-conservation of the in-plane momentum in the elec- 
tron excitation process. Hence, the phototransitions 
should include the participation of the "third body". In 
our case the impurities play the role of this agent. The 
excitation probability including the impurity scattering 
is determined by the second-order transition amplitude. 
The needed term arises from the interference of ampli- 
tudes caused by the Ez and in-plane components of the 
electric field. The draft of the transitions is depicted in 
Fig.l. 

In the second order of the interaction, the transition 
probability is 



w- 



ph 



n,p;n ,p 




V tJI^ 

^ n,p;ni,p' <~J ni,p\n' ,p' 



"^n,p;ni,p ''ni,p:n',p' 



+ uj) 77 + i(e„i^„ - cj) 



^^,p;ni ,p' ,p;n',p' 



^n,p;ni ,p^^i ,p;n' ,p' 

V + i{Sm,n + w) 



^n' ,p' 



K^n.p 

Here Em, 

erage over impurities configuration 



u); (77 = +0). 

~ e„; angular brackets denote the av- 



Using relations 



n,p;n',p' {U n' .p' :n ,p) : 14i,p:n',p' — (14i' ,p';n,p) it iS 
ph 



n ,p ;n,p' 



easy to prove that W^_p.„,_p 

The denominators in Eq. ^ have their resonance with 
the field frequency independently from the electron mo- 
mentum. At the same time, the resonance in the final 
state is absent due to non-conservation of the in-plane 
momentum. 

Eq.® can be rewritten in the form (E = (En, E'z)): 




Vp'E||(5„j „' 



,E* 

ni.n' z 



V'' E* 



Tj -\- i (^m ,n' 
\ 2 



V + i{enun - W) 



Vp'E||(5„i^„/ 



K.mEz 



ni ,p;n' ,p' 



K^^n'Ez 
\ 2 



Wr 



ni ,p;n' ,p' 



It is evident that the contribution to photogalvanic ef- 
fect is given by not the total transition probability W^^^ 
but only its odd in p, p' part. For this part, we have the 
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following expression: 



ph 



n.p-.n' ,p' 
2 



Re 



Er-p;"%p'(p'-P)Ei 



*^n,p;ni,p'^n',ni^ 



E,V* , , 



'5(en,p - En'.p' + 1^) + (Re 



V* /v"", E* 



K,P;«',P'(P-P')E|| 



^ni^n-^z •^ni,p;n',p' 



l5(en,p - en',p' - w) >. 



(10) 



Kinetic equation Eq.(j4]) can be transformed to 

f (1) _ r 

rJ-n,p ~ ^n,pi 



(11) 



where T„(p) is the intra-subband transport relaxation 
time and t„ _„(p) is the time of transition from the state 
{n, p) to all states of the subband (— n). These values are 
determined by 



^(.P) 



P' 



Vr 



n,p;n,p' | /^(^n.p ^n,p')(^ 



^2 



^ (^n,p ^— n,p' ) 



''n,p; — n,p'| ^^(^n,p ^— n.p') ^2 



12) 



Solving Eg. (1111) we find (argument p is omitted): 



(13) 



The expressions for r„(p), t„ _„(p) and W,^_p.„/_p/ con- 
tain correlators of the form (Vn^p-n' ,p'Vm,p;m' ,p''^ ■ In the 

case of impurities situated in layer z = zq (r^ — (pj, zq)) 
the function F(r) reads 



Vir) = ^ wqe-«l^-^»l exp (-*q(p - ft)), 



(14) 



where Uq is the 2D Fourier component of the impurity 
center potential. For example, for unscreened Coulomb 
center Uq = 2i:e'^ j nqS {k is the background dielectric 
constant). Correlators are given by 



K,p;n',p'K;,p;™',p') = n,5/dzdz>p_pf e-^(2-0-.-.') 



Here rig is the areal density of scatterers. We suppose 
that the electron wavelength is larger than d. In this 
approximation one can find from Eq. (|15p : 



^n,p;n' ,p'yra,p:m' ,p'/ ~ ^s'S'lUp^p'j 6 " X 



Matrix elements z„„/ should be estimated for specific 
wave functions. For simplicity, we will use the wave func- 
tions of two delta-functional wells in the tight-binding 
approximation. The seed states with energies Eq ± A.q/2 
can be written as 



Xi,2 = 



(17) 



In basis m X+ = (l,/3)/yT+^, X- = 
(/3, —1)/ •\/l + /3^, where /? is the mixing amplitude. The 
corresponding states energies are e± = Eq ± A/2, A = 
a/Aq -I- 4iQ, where to ~ £oe~'"' is a hopping amplitude 
between wells. For quantity /3 we have /3 = 2to/(A + Aq). 

The matrix elements of z are z++ = — z = d(l — 

/?2)/(2(l-h/32)), z+_=/3d/(l-H/32). 

Inserting Eq.(fT6|) in Eq. ([T2l) we get the expressions for 
T+ w T_ = r and a small difference 1/t_ — 1/t+: 



— = mn^ 

T 



Zn 



%2 + 2pq)^ 



= m(z++ — z )ns X 

r_ T+ 



TT 



%2 + 2pq)%, 



(18) 



where Uq = Suq. From Ea. ((T5|) it is seen that r„^_„ ^ r„ 
and, so, Eq.([T3]) can be simplified 



f(i) - G T 

Jn ~ In- 



(19) 



Further we will consider the resonance situation when 
frequency w is close to A. Smallness A, as compared 
to the Fermi energy ep — pj?/2m [pp being the Fermi 
momentum) leads to approximate expressions for Gj^ ~ 
-G_, 



nse'^Az1_ 



(0) 



TTW^ (A — + "ff 



J dq(q • Rc{E^^E:))q\uq\'e-^'^'°S{q^ + 2pq). (20) 

As a result, for the current of photogalvanic effect, we 
have 



J = 



An^muj^ (A - uj)^t^ + 1 



Re(E||£;;) X 



dp 



9/i 



dq\u. 



~ \2-2qzo^3 



q^<5(g^ + 2qp). (21) 



Eq. (|2ip has a resonant character with the resonance at 
^Xn{z)xn'{z)Xm{z')Xrn'{z'). (15) w — A. This rcsonancc results from the intermediate 
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state for transition due to the parallelism (equidistance) 
of subbands. The resonance is smeared due to scatter- 
ing, e.g., by impurities. To include this smearing, the 
infinitesimal rj was replaced by finite relaxation rate 1 /r 
which can be estimated from mobility. This leads to the 
flniteness of the current at the point of resonance. 

At temperature T = the latter expression is simpli- 
fied, and we obtain the final result for the required value: 



4e^(z++ - z )Az^„eFT 

7ra;2rf2((A-a;)2T2 + l) 



Rei^lE*)F, (22) 



where we introduced a dimensionless quantity F 



2pF 



1 



(23) 



In the specific case oippZQ ^ 1 Eq. 



is reduced to 



(24) 



If the scattering is determined by the charged non- 
screened impurities F = d'^ /Az^. 

In the model of two (5-like wells Eq. gives For linear 
polarized wave 



4e3(i/3(l - I3^)n 



i(l-t-/32)2A (A/fi-a;)2T2 + l 



sin (261) F, 



(25) 

where Eq is the amplitude of the electric field, 6 is the 
angle between the field and the normal to the system. 

It should be emphasized that the current contains the 
linear response only. This distinguishes the quantum 
double-well result from the simple classical model of the 
effect considered in the previous section. 

Let us estimate the value of the effect. Considering /3 
as a free parameter we can choose (i = \/2— 1 to maximize 
the /3-dependent factor in Eq.(l25]) /3(1 - ^2)/(l -h /3^)2 = 



1/4. The optimum for PGE observation corresponds to 
Lu — Is. and 9 = 7r/4. Choosing the typical values for 
GaAs/AlGaAs double quantum wells d — h ■ 10~^ cm, 
zo = S-lO^^cm, 6^=20 meV (n = 6.2- 10"cm-2), A=0.1 
meV, T = 4 • 10^^^ s, £^o = 1 V/cm we find for this opti- 
mal situation j « 3.6 /iA/cm that is a quite measurable 
value. 

It should be emphasized that the initial and final states 
in the transition can belong to the different or the same 
subbands. The resonant behavior results from the reso- 
nance on the intermediate state rather than the energy 
conservation in the final states, because the conservation 
law for the phototransition with the participation of im- 
purity scattering does not give a fixed frequency for the 
transition. The sharpness of the resonance is conditioned 
by the equidistance of the energy bands in a 2D well. 
IV. CONCLUSIONS 

We found the stationary current along a double-well 
system affected by the linear-polarized far-infrared wave. 
The stationary current originates from the periodic vi- 
bration of electrons between two non-equivalent quan- 
tum wells caused by the normal component of the al- 
ternating electric field with synchronic in-plane accelera- 
tion/deceleration by the in-plane component of the elec- 
tric field. The linear photogalvanic effect needs vertical 
asymmetry of the quantum well. The effect has the peak 
resonant structure connected with the parallel subbands 
of the double quantum well. The resonant frequency can 
be easily tuned by the application of the gate voltage. 
The optimal range of frequencies is lO^^H- 10^'^ s^^. 
The predicted value of the current is experimentally mea- 
surable. 
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